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analysis'"*  yields  the  bounds  >  K  >  Ymm  where,  to  first  order  in 
the  small  S, 


Kmax.min 
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Again,  these  quantities  will  be  more  sensitive  if  the  moments  of 
inertia  are  close  than  if  they  are  widely  spaced.  For  a  vehicle  with 
nominal  inertias  satisfying  Ii  =  2Ij  and  l2  =  1.5 13,  and  the  1  %  nor¬ 
malized  inertia  perturbations  considered  earlier,  arctan  y  can  vary 
from  31.8  to  38.7  deg,  a  total  range  of  6.9  deg.  Using  Eqs.  (13)  and 
(14)  together  allows  the  possible  variations  in  absolute  separatrix 
directions  to  be  quantified  as  a  function  of  the  inertia  error  bound  S, 
thus  providing  the  desired  insight  into  sizing  Sy  as  a  function  of  S. 

Conclusions 

This  Note  has  derived  a  design  procedure  that  prevents  a  space¬ 
craft  with  a  stuck-on  thruster  from  experiencing  a  large  net  linear 
acceleration  and,  hence,  a  significant  perturbation  to  its  orbit.  It  was 
shown  that  such  an  acceleration  is  most  likely  to  arise  for  a  thruster 
that  generates  a  torque  that  is  nominally  about  the  intermediate  axis 
of  an  asymmetric  vehicle,  but  that  is  in  reality  slightly  offset  from 
it.  The  simple  technique  that  was  developed  to  avoid  this  makes  use 
of  a  small  shift  in  thruster  position.  Expressions  were  derived  to 
allow  this  shift  to  be  sized  to  produce  robust  results  despite  small 
uncertainties  in  the  geometry  and  mass  properties  of  the  vehicle. 

Acknowledgments 

This  work  was  supportedinpartby  the  Automation,  Robotics,  and 
Simulation  Division  at  NASA  Johnson  Space  Center  under  Grant 
NAG9-874,  with  Technical  Monitor  Keith  Grimm.  The  authors  also 
wish  to  thank  the  Associate  Editor  and  reviewer  for  many  helpful 
suggestions. 

References 

^Gemini  Summary  Conference,  NASA  SP-138,  Feb.  1967,  pp.  52,  165. 

^Hope,  A.  S.,  and  Middour,  J.  W.,  “Clementine  Mission  Maneuver  Per¬ 
formance  and  Navigation  Control  for  the  Spinning  Phase,”  Guidance  and 
Control  1995,  Advances  in  the  Astronautical  Sciences,  Vol.  88,  Univelt,  San 
Diego,  1995,  pp.  463-476. 

^WIRE  Mishap  Investigation  Board  Report,  NASA,  June  1999,  pp.  14, 
15. 

"*Leimanis,  E.,  The  General  Problem  of  the  Motion  of  Coupled  Rigid 
Bodies  About  a  Fixed  Point,  Springer- Verlag,  Berlin,  1965,  pp.  136-178. 

^Chobotov,  V.  A.,  Spacecraft  Attitude  Dynamics  and  Control,  Krieger, 
Malabar,  FL,  1991,  pp.  26-29. 

^Randall,  L.  A.,  Longuski,J.  M.,  and  Beck,  R.  A.,  “Complex  Analytic  So¬ 
lutions  for  a  Spinning  Rigid  Body  Subject  to  Constant  Transverse  Torques,” 
Astrodynamics  1995,  Advances  in  the  Astronautical  Sciences,  Vol.  90,  Pt.  I, 
Univelt,  San  Diego,  1995,  pp.  593-607. 

^Tsiotras,  P,  and  Longuski,  J.  M.,  “A  Complex  Analytic  Solution  for 
the  Attitude  Motion  of  a  Near-Symmetric  Rigid  Body  Under  Body-Fixed 
Torques,”  Celestial  Mechanics  and  Dynamical  Astronomy,  Vol.  51,  No.  3, 
1991,  pp.  281-301. 

^Williams,  T.  W.,  and  Tanygin,  S.,  “Dynamics  of  a  Near-Symmetrical 
Spacecraft  Driven  by  a  Constant  Thrust,”  Proceedings  of  6th  AiAS/AIAcA 
Spaceflight  Mechanics  Meeting,  Advances  in  the  Astronautical  Sciences, 
Vol.  93,  Pt.  II,  Univelt,  San  Diego,  1996,  pp.  925-944. 

^Fivneh,  R.,  and  Wie,  B.,  “New  Results  for  an  Asymmetric  Rigid  Body 
with  Constant  Body-Fixed  Torques,”  Journal  of  Guidance,  Control,  and 
Dynamics,  Vol.  20,  No.  5,  1997,  pp.  873-881. 

'"williams,  T.  W.,  and  Tanygin,  S.,  “Dynamics  of  an  Asymmetric  Space¬ 
craft  Driven  by  a  Constant  Thrust:  Perturbation  Studies,”  Proceedings  of  7th 
AAS/AIAA  Spaceflight  Mechanics  Meeting,  Advances  in  the  Astronautical 
Sciences,  Vol.  95,  Pt.  I,  Univelt,  San  Diego,  CA,  1997,  pp.  211-229. 

"Bracewell,  B.  N.,  and  Garriott,  O.  K.,  “Rotation  of  Artificial  Earth 
Satellites,”  Vatwre,  Vol.  182,  No.  4638, 1958,  pp.  760-762. 

'^Kaplan,  M.  H.,  Modem  Spacecraft  Dynamics  and  Control,  Wiley,  New 
York,  1976,pp.60,61. 

'"Tanygin,  S.,  “The  General  Problem  of  Attitude  Motion,”  Ph.D. 
Dissertation,  Dept,  of  Aerospace  Engineering,  Univ.  of  Cincinnati, 
Cincinnati,  OH,  Feb.  1999,  pp.  14-19. 

'"*Wilkinson,  J.  H.,  The  Algebraic  Eigenvalue  Problem,  Oxford  Univ. 
Press,  Oxford,  England,  U.K.,  1965,  pp.  69-75. 


G.  Song* 

University  of  Akron,  Akron,  Ohio  44325-3903 
and 

S.  P.  Schmidt^  and  B.  N.  Agrawal*^ 

Naval  Postgraduate  School,  Monterey,  California  93943 

I.  Introduction 

OSITIVE  position  feedback  (PPE)  control  was  introduced  by 
Goh  and  Caughey  in  1985  to  control  vibrations  of  large  flex¬ 
ible  space  structures.'  A  PPE  controller  has  several  distinguished 
advantages  as  compared  to  then  widely  used  velocity  feedbackcon- 
trol  laws.  It  is  insensitive  to  spillover"  where  contributions  from 
unmodeled  modes  affect  the  control  of  the  modes  of  interest."  "'  As 
a  second-order  low-pass  filter,  a  PPE  controller  rolls  off  quickly  at 
high  frequencies  and  is  well  suited  to  controlling  the  lower  modes 
of  a  structure  with  well  separated  modes."  "  Because  of  these  ad¬ 
vantages,  PPE  controller  along  with  smart  materials,  in  particular 
PZT  (lead  zirconate  titanate)  type  of  piezoelectric  material,  has  been 
applied  to  many  flexible  systems  to  achieve  active  damping."""  The 
design  of  a  PPE  controllerrequires  the  natural  frequency  of  a  struc¬ 
ture.  In  practice,  the  structural  natural  frequency  may  not  be  known 
exactly  or  it  may  vary  with  time.  When  the  frequency  used  in  the  PPE 
controller  is  different  from  that  of  the  structure,  the  performance  of 
the  PPE  control  will  adversely  affected.  Despite  that  PPE  control 
is  widely  researched  in  literature,  robustness  study  of  PPE  control 
when  natural  frequency  is  inexactly  known  is  not  reported.  This 
motivates  the  authors  to  conduct  experimental  study  of  robustness 
of  PPE  control  in  active  vibration  suppression  of  a  smart  flexible 
structure. 

II.  PPF  Control 

PPE  control  requires  that  the  sensor  is  collocated  or  nearly  collo¬ 
cated  with  the  actuator.  In  PPF  control  structural  position  informa¬ 
tion  is  fed  to  a  compensator.  The  output  of  the  compensator,  mag¬ 
nified  by  a  gain,  is  fed  directly  back  to  the  structure.  The  equations 
describing  PPF  operation  are  given  as 

^  -E  2^0}^  -\-  =  Gofrj  (1) 

rj  +  2l;cCo]i]  +  co]i!)  =  co]^  (2) 

where  |  is  a  modal  coordinate  describing  displacementof  the  struc¬ 
ture,  f  is  the  damping  ratio  of  the  structure,  w  is  the  natural  frequency 
of  the  structure,  G  is  a  feedback  gain,  rj  is  the  compensator  coordi¬ 
nate,  fc  is  the  compensator  damping  ratio,  and  coc  is  the  frequency 
of  the  compensator.  For  the  closed  loop  being  stable,  0  <  G  <  1 
(Ref.  4). 

To  illustrate  the  operation  of  a  PPF  controller,  assume  a  single 
degree-of-freedom  vibration  of  the  beam  in  the  form  of 

m  =  (3) 

the  output  of  the  compensator  at  the  steady  state,  provided  the 
closed-loop  system  is  stable,  will  be 

rj(t)  =  (4) 

Received  2  April  200 1 ;  revision  received  30  July  2001 ;  accepted  for  pub¬ 
lication  7  September  2001.  This  material  is  declared  a  work  of  the  U.S. 
Government  and  is  not  subject  to  copyright  protection  in  the  United  States. 
Copies  of  this  paper  may  be  made  for  personal  or  internal  use,  on  condi¬ 
tion  that  the  copier  pay  the  $  1 0 .00  per-copy  fee  to  the  Copyright  Clearance 
Center,  Inc.,  222  Rosewood  Drive,  Danvers,  MA  01923;  include  the  code 
0731-5090432  $10.00  in  correspondence  with  the  CCC. 

*  Assistant  Professor,  Department  of  Mechanical  Engineering. 

'  Student,  Department  of  Aeronautics  and  Astronautics. 

^Professor,  Department  of  Aeronautics  and  Astronautics.  Associate  Fel¬ 
low  AIAA. 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

2QQ  j  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2001  to  00-00-2001 

4.  TITLE  AND  SUBTITLE 

Experimental  Robustness  Study  of  Positive  Position  Feedback  Control  for 
Active  Vibration  Suppression 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School, Department  of  Aeronautics  and 

Astronautics, Monterey, CA, 93943 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

4 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


180 


J.  GUIDANCE,  VOL.  25,  NO.  1 :  ENGINEERING  NOTES 


In  Eq.  (4),  the  magnitude  /J  is  given  as 

A(ft)/wc) 

P  /  2 - 2 

y  (l  -  wlY  +  [2fc(wM)]^ 


where  A  =  a{a)clco). 

In  Eq.  (4),  the  phase  angle  4>  is  given  as 


4>  =  tan  ' 


2fc(a>/uic) 

1  -  (cu^/w?) 


(5) 


When  the  structural  modal  frequency  is  much  lower  than  the  com¬ 
pensator  natural  frequency,  the  phase  angle  (p  approaches  zero  ac¬ 
cording  to  Eq.  (5).  Substituting  Eq.  (4)  with  0  =  0  into  Eq.  (1) 
results  in 


$  +  2f  cui  -b  (w"  -  =  0  (6) 

It  is  clear  from  Eq.  (6)  that  the  PPE  compensatorin  this  case  results 
in  the  stiffness  term  being  decreased,  which  is  called  active  flexibil¬ 
ity.  When  the  compensator  and  the  structure  have  the  same  natural 
frequency,  it  can  be  derived  from  Eq.  (6)  that  the  phase  angle  0 
approaches  7r/2.  Substituting  Eq.  (4)  with  0  =  7r/2  into  Eq.  (1),  the 
structural  equation  becomes 

|-b(2fu;-bG/Scu)i-bu;"|  =0  (7) 

Equation  (7 )  shows  that  the  PPE  compensator  in  this  case  results  in 
an  increase  in  the  damping  term,  which  is  called  active  damping. 
When  the  structure  frequency  is  much  greater  than  that  of  the  com¬ 
pensator,  the  phase  angle  0  approaches tt.  SubstitutingEq.  (4)  with 
0  =  jr  into  Eq.  (1)  results  in 

I  -b  2f  -b  {ofl  -b  =  0  (8) 


III.  Experimental  Setup 

A  schematic  of  the  equipment  setup  for  vibration  control  using 
PPE  is  shown  in  Eig.  2.  A  cantileveredaluminum beam  (its  properties 
are  shown  in  Table  1 )  is  used  as  the  object  for  vibration  control.  The 
beam  has  a  PZT  sensor  and  three  PZT  actuators.  The  PZT  sensor 
actually  consists  of  a  pair  of  PZT  patches,  one  on  each  side  of  the 
beam.  Each  PZT  actuatoralsoconsistsof  apair  of  PZT  patches,  one 
on  each  side  of  the  beam.  The  PZT  sensor  can  be  considered  nearly 
collocated  with  PZT  actuator  1.  Properties  of  the  PZT  are  shown  in 
Table  2.  The  aluminum  beam  is  clamped  such  that  its  length  was 


Table  1  Cantilevered  beam  properties 
(aluminum  beam  type:  7075  T-6) 


Quantity 

Description 

Value 

tb 

Beam  thickness 

1.95  mm 

Wb 

Beam  width 

76.2  mm 

L 

Beam  length 

1.176m 

Pb 

Beam  density 

2800  kg/m^ 

Eb 

Young’s  modulus 

7.1  X  10'°N/m^ 

Table  2  Properties  of  the  piezoceramic  actuators 
and  sensor  [Type:  PZT-5A  (Navy  type  II)] 


Quantity 

Description 

Value 

d-ii 

Lateral  strain  coefficient 

1.8  X  10““  Coul/N 

Ep 

Young’s  modulus 

6.3  X  10“  N/m^ 

el 

Absolute  permittivity 

1.5  X  10“*  Farad/m 

Pp 

PZT  density 

7.7  X  10^  kg/m^ 

h 

PZT  actuator  and  sensor  thickness 

0.5  mm 

Wp 

PZT  actuator  and  sensor  width 

38.1  mm 

La 

Length  of  PZT  actuators 

63.5  mm 

L, 

Length  of  PZT  sensor 

25.4  mm 

It  is  clear  that  from  Eq.  (8)  that  the  PPE  compensatorin  this  case 
results  in  an  increase  in  the  stiffness  term,  which  is  called  active  stiff¬ 
ness.  To  achieve  maximum  damping,  01^  should  be  closely  matched 
to  (d.  Also,  any  structural  natural  mode  below  01^  will  experience 
increased  flexibility. 

The  effect  of  the  damping  ratio  is  discussed  as  follows.  Larger 
values  of  the  damping  ratio  f  ^  will  result  in  a  les  s  steep  slope  thereby 
increasing  the  region  of  active  damping.  Eigure  1  shows  the  bode 
plot  for  fc  =  0.5  and  for  =  0. 1 .  The  difference  in  the  slopes  of  the 
phase  angle  can  be  easily  seen.  A  larger  value  of  ensures  a  larger 
region  of  active  damping  and  therefore  will  increase  the  robustness 
of  the  compensator  with  respect  to  uncertain  modal  frequencies. 
However,  it  is  expected  to  result  in  slightly  less  effective  damping 
and  result  in  increased  flexibility  at  lower  modes  as  a  tradeoff. 


PZT 

Sensor 

PZT  PZT  PZT 

Actuator#]  Actuator#2  Actuator#3 


Fig.  2  Experimental  setup  schematic. 


Gain  dB 


Gain  dB 


Fig.  1  Bode  plot  for  PPF  with  C^c  =  0.5  (left)  and  C^c  =  0.1  (right). 


J.  GUIDANCE,  VOL.  25,  NO.  1 :  ENGINEERING  NOTES 


181 


parallel  to  the  granite  table  below  it.  This  allowed  the  bending  to 
be  strictly  in  the  horizontal  plane.  A  digital  control  system,  called 
Modular  Control  Patch,  is  used  to  implement  vibration  suppression 
algorithms.  It  uses  a  TMS320C30  microprocessor.  A  digital  data 
acquisition  system  using  a  TMS320C40  digital  signal  processor 
was  used  to  record  the  experimental  data. 

IV.  Experimental  Procedure 

Both  open-  and  closed-loop  tests  were  performed.  All  tests  were 
started  by  manually  exciting  the  beam.  This  was  a  simple  and  an 
effective  method  to  excite  the  beam.  For  single-mode  vibration  sup¬ 
pression  tests  can  be  run  with  eitherall  three  actuatoroperationafor 
only  the  first  actuator  operational.  For  multimode  suppression  only 
the  first  PZT  actuator  was  used.  Because  of  modal  shapes  for  the 
higher  modes  of  the  flexible  beam,  the  second  and  the  third  actuators 
adversely  impact  damping  of  higher  modes  when  the  feedback  strain 
information  from  the  only  PZT  sensor  is  used.  It  is  clear  from  Fig.  2 
that  the  PZT  sensor  is  not  collocated  or  near  collocated  with  these 
two  actuators.  For  each  test  data  were  recorded  for  a  time  interval 
of  15  s  after  beam  excitation.  This  allowed  ample  time  to  measure 
damping  effects.  The  experimentaldata  were  then  processed  to  show 
effectiveness  of  the  tested  control  algorithm.  A  fast  Fourier  trans¬ 
form  was  performed  to  provide  a  power  spectral  density  (PSD)  plot 
of  the  beam  response.  PSD  gives  a  measure  of  signal  energy  level  at 
different  frequencies.  A  comparison  of  the  ratios  of  the  last-second 
modal  energy  level  in  decibels  to  the  initial  one  provides  an  indi¬ 
cation  of  the  damping  effectiveness  on  this  particular  mode.  Also, 
a  direct  comparison  of  the  modal  energy  level  drop  with  that  of 
an  open-loop  response  can  indicate  the  effectiveness  of  the  control 
algorithm. 

Figure  3  shows  the  PSD  plots  for  a  multimode  open-loop  vibra¬ 
tion.  The  solid  line  is  for  the  first  second  of  the  15-s  test  and  the 
dashed  line  for  the  last  second.  A  program  was  written  to  iden¬ 
tify  the  modes  excited  and  to  compute  the  difference  between  the 
initial  and  final  energy  level  in  the  decibel  unit  at  the  identified 
modal  frequencies.  The  following  data  show  energy  level  drops  in 
decibels  for  the  first  four  modes  in  the  15-sfree  vibration:  first  mode 
(1.33  Hz),  9.52  dB;  second  mode  (7.1  Hz),  22.38  dB;  third  mode 
(19.0  Hz),  48.98  dB;  fourth  mode  (38.2  Hz),  61.94  dB.  It  is  clear 
that  vibrations  of  the  third  and  fourth  modes  quickly  damp  out. 
The  first  and  second  modes  become  the  major  concern  for  vibration 
suppression. 

V.  Experimental  Robustness  Study  of  PPE 
in  Single  Mode  Vibration  Suppression 

One  of  the  documented  drawbacks  of  PPF  is  that  its  design  re¬ 
quires  the  knowledge  of  modal  frequency.  If  the  targeted  frequency 
isalteredoris  simply  miscalculated,  PPF  damping  would  be  severely 


Fig.  4  Robustness  results  for  different  (^c- 


affected.  To  test  the  robustness  of  PPF  with  respect  to  the  modal  fre¬ 
quency,  experiments  of  PPF  with  inaccurate  modal  frequency  from 
25  to  400  %  of  the  nominal  value  were  conducted. Experiments  with 
different  value  of  were  also  tested.  Three  values  were  used  for 
fc!  1,  0.5,  and  0.25.  This  set  of  experiments  was  carried  out  to  study 
the  impact  of  compensator  damping  ratio  on  compensator  robust¬ 
ness  and  on  possible  increased  flexibility  at  lower  modes.  From  the 
analysis  in  Sec.  II,  increasing  the  value  of  the  damping  ratio  will 
provide  a  wider  frequency  range  for  active  damping,  but  it  is  ex¬ 
pected  that  its  effectiveness  at  the  target  frequency  will  be  reduced. 
In  these  experiments  only  PZT  actuator  one  is  used  along  with  the 
PZT  sensor.  The  PPF  controller  employs  G  =  0.086. 

The  experimental  results  of  suppression  of  first  modal  vibration 
using  the  PPF  controlareplottedinFig.4.  In  this  figure  the  percent¬ 
age  values  indicate  percentage  increase  of  the  modal  energy  drop  in 
decibels  with  PPF  control  as  compared  to  that  in  the  case  of  free  vi¬ 
bration.  It  can  be  seen  from  this  figure  that  the  robustness  increases 
as  fc  is  increased.  This  is  as  expected.  The  effectiveness  of  the  PPF 
at  the  target  frequency  of  1 .3  Hz  was  only  slightly  reduced  when 
changing  damping  ratio  from  0.25  to  1.0.  This  observation  recom¬ 
mends  using  higher  value  of  in  the  PPF  compensator  design  to 
achieve  robustness  because  effectiveness  at  the  targeted  mode  will 
not  be  much  affected  by  a  higher  PPF  control  showed  positive 
damping  on  all  frequencies  tested  but  dropped  off  rapidly  when  go¬ 
ing  below  75%  of  the  targeted  frequency  or  200%  above  it.  Overall, 
PPF  control  is  robust  to  modal  frequency  variations.  In  addition, 
increased  flexibility  was  not  measurable  even  at  300%  above  the 
fundamental  frequency.  An  optimal  for  a  given  structure  would 
depend  on  how  accurate  the  modes  are  known  and  how  much  they 
wouldbe  expected  to  change,  but  a  fc  of  at  least  0.5  shouldbe  chosen 
for  robustness. 

VI.  Experimental  Robustness  Study  of  PPE 
in  Multimode  Vibration  Suppression 

PPF  controls  were  tested  for  multimode  vibration  suppression. 
Two  PPF  filters  are  used.  The  first  PPF  filter  targeted  the  first  mode, 
and  the  second  PPF  filter  targeted  the  second  mode.  The  two  PPF 
filters  are  in  a  parallel  connection.  Only  actuator  1  along  with  the 
senor  is  used.  Experiments  were  conducted  to  investigate  the  ro¬ 
bustness  of  this  type  of  controller.  Table  3  shows  the  test  results. 
In  this  table  the  percentage  values  indicate  percentage  increase  of 
the  modal  energy  drop  in  decibels  with  PPF  control  as  compared  to 
that  in  the  case  of  free  vibration,  =  0.5  is  used  for  both  filters.  For 
the  first  filter  G  =  0.086,  and  for  the  second  filter  G  =  0.003.  Both 
compensator  frequencies  are  moved  progressively  higher  than  the 
targeted  natural  frequency.  As  seen  from  the  table,  the  first  mode 
damping  falls  off  much  quicker  than  the  second  mode.  Part  of  the 
reason  for  this  is  the  increased  stiffness  region  of  the  first  filter  is 
moving  closer  to  the  second  mode,  thereby  helping  to  increase  the 
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Table  3  Robustness  results  for  two  PPF 
filter  combination 


PPF  modal  frequency, 

Hz 

First  mode 
decibel  drop 
(%  increase) 

Second  mode 
decibel  drop 
(%  increase) 

(jOci  =  1.3  X  1.0=  1.3 

0)^2  =  7.1  X  1. 0  =  7.1 

58.33(513%) 

44.00  (97%) 

cuci  =  1-3  X  1.25  =  1.625 
Mc2  =  V  l  X  1.25  =  8.875 

42.16(343%) 

38.10(70%) 

Mci  =  1.3  X  1.5  =  1.95 
cuc2  =  7.1  X  1.5  =  10.65 

36.14(280%) 

36.35  (62%) 

cuci  =  1-3  X  1.75  =  2.275 
cUc2  =  7.1  X  1.75  =  12.45 

30.11  (216%) 

32.99  (47%) 

coci  =  1.3  X  2  =  2.6 

C0c2  =  l-l  X  2=  14.2 

17.88  (88%) 

26.57(19%) 

dB 


Fig.  5  PSD  plot  for  two  PPF  filters  (compensator  frequencies  are  1.5 
times  the  targeted  modal  frequencies.). 


damping  effect.  Just  as  in  the  single  PPF  case,  the  PPF  combination 
shows  good  robustness  for  both  modes.  These  results  suggest  that 
along  with  being  robust  the  two  PPF  filter  combination  is  effective  in 
damping  multiple  modes  over  a  range  of  frequencies .  A  PSD  plot  for 
the  same  controllergraphically  illustrating  damping  effectiveness  is 
shown  in  Fig.  5. 

VII.  Conclusions 

This  research  presents  the  experimental  results  robustness  study 
of  vibration  suppression  of  a  flexible  structure  using  PPF  control. 
The  flexible  structure  is  a  cantilevered  beam  with  PZT  sensors  and 
PZT  actuators.  PPF  controls  were  implemented  for  single-mode 
vibration  suppression  and  for  multimode  vibration  suppression.  Ex¬ 
periments  found  that  PPF  control  is  robust  to  frequency  variations 
for  single-mode  and  for  multimode  vibration  suppressions. 
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Approximate  Equations  for 
the  Coplanar  Restricted 
Three-Body  Problem 
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Introduction 

HE  fuel-optimal  rendezvous  problem  for  a  spacecraft  with  a 
satellite  in  a  central  force  field  has  been  the  subject  of  several 
papers  in  the  last  decade.  In  this  Note  we  address  the  rendezvous 
problem  of  a  spacecraft  with  a  moon,  comet,  or  asteroid  (whose 
gravitational  field  is  not  negligible)  in  the  central  force  field  of  a 
(large)  third  body.  In  this  setting  we  derive  reduced  equations  for 
the  motion  of  the  spacecraft  in  the  vicinity  of  the  smaller  body 
(moon)  in  three  dimensions.  We  present  a  generalization,  under 
some  restrictions,  to  the  Jacobi  integral.^ 

Clearly  the  problem  we  consider  here  is  closely  related  to  the 
restricted  three-body  problem  in  two  dimensions,^-*  which  deals 
with  the  motion  of  a  body  of  negligible  mass  in  the  gravitational 
field  of  two  other  celestial  bodies.  However  in  this  setting  all  three 
bodies  are  assumed  to  remain  always  in  one  plane. 

To  put  our  results  in  proper  perspective,  we  note  that  Edelbaum^ 
made  one  of  the  first  attempts  to  address  the  rendezvous  problem 
between  a  spacecraftand  a  satellitein  a  near  circularorbit.  A  simpler 
model  was  derived  by  Clohessy  and  Wiltshire,'”  whose  equations 
can  be  found  now  in  books  on  orbital  mechanics.''  Carter  and  Humi,' 
Humi,^  and  Carter  and  Brient*  found  some  analytical  solutions  for 
the  rendezvous  problem  when  the  satellite  is  in  a  general  Keplerian 
orbit  and  derived  a  generalizatior?  of  these  equations  in  the  pres¬ 
ence  of  a  general  central  force  field.  Currently  an  effort  is  under 
way  to  include  the  effect  of  drag  (linear  and  quadratic)  in  the  ren¬ 
dezvous  equations.”  (For  a  more  extensive  list  of  contributions  to 
the  rendezvous  problem,  see  the  reference  lists  in  Refs.  1  and  6). 

Derivation  of  the  Reduced  Equations 

In  an  inertial  coordinate  system  whose  origin  is  at  the  center  of 
the  central  body  E,  we  let  R,  p  denote  the  position  of  the  moon 
and  the  spacecraft,  respectively,  and  r  the  relative  position  of  the 
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